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We express Einstein's field equations for a spherically symmetric ball of general fluid such that they 
are conducive to an initial value problem. We show how the equations reduce to the Vaidya spacetime 
in a non-null coordinate frame, simply by designating specific equations of state. Furthermore, 
this reduces to the Schwarzschild spacetime when all matter variables vanish. We then describe 
the formulation of an initial value problem, whereby a general fluid ball with vacuum exterior is 
established on an initial spacelike slice. As the system evolves, the fluid ball collapses and emanates 
null radiation such that a region of Vaidya spacetime develops. Therefore, on any subsequent 
spacelike slice there exists three regions; general fluid, Vaidya and Schwarzschild, all expressed 
in a single coordinate patch with two free-boundaries determined by the equations. This implies 
complicated matching schemes are not required at the interfaces between the regions, instead, one 
simply requires the matter variables tend to the appropriate equations of state. We also show the 
reduction of the system of equations to the static cases, and show staticity necessarily implies zero 
"heat flux" . Furthermore, the static equations include a generalization of the Tolman-Oppenheimer- 
Volkoff equations for hydrostatic equilibrium to include anisotropic stresses in general coordinates. 

PACS numbers: 04.20.-q, 04.40.Nr, 04.20.Jb 



I. INTRODUCTION 



Full analytic studies of gravitational collapse are intrin- 
sically difficult due to the non-linearities in the Einstein 
field equations (EFEs). Indeed, in the study of gravita- 
tional collapse, a standard approach is to approximate 
the region exterior to the collapse to be vacuum (see 
for example [J, Q and references therein), which under 
the condition of spherical symmetry is the Schwarzschild 
spacetime. While this simplification reduces the math- 
ematics to tangible systems, it is physically insufficient. 
The next approximation beyond this is to allow for some 
form of null radiation to emanate from the collapsing re- 
gion [1, 0, H, @, 13 • I n this way, the number of baryons 
are conserved within the collapsing region, however en- 
ergy may still be dissipated. In spherical symmetry, a re- 
gion of outgoing null radiation is described by the Vaidya 
spacetime. 

The Vaidya spacetime is a solution of EFEs that con- 
tains incoherent null radiation; non- interacting particles 
moving along null geodesies. References d, H, [l(| have 
shown the Vaidya spacetime to be an extension of the 
dust filled Lcmaitrc-Tolman-Bondi (LTB) spacetime. To 
do this, they proposed a LTB interior region, and at some 
finite radius the timelike coordinate axis tilted over to be- 
come null. While the mathematics of this method works 
well, the physical realization it brings is intuitively awk- 
ward. This is due to the exterior region of the spacetime 
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still being in a null coordinate frame, implying the set-up 
is not conducive to an initial value problem. 

To overcome this we recast the entire problem in a 
non-null coordinate frame such that it can be used as an 
initial value problem. In this frame, we show the incoher- 
ent null radiation of the Vaidya spacetime takes the form 
of a general fluid with very specific equations of state. 
Therefore, rather than our interior region being dust, we 
formulate the interior region as a general spherically sym- 
metric fluid such that the energy momentum terms are 
able to smoothly match to a Vaidya exterior. Therefore, 
this article generalizes our earlier work [Til . Il2| by formu- 
lating a completely general, spherically symmetric fluid 
utilizing the 3 + 1 formalism for general relativity which 
thus makes it suitable for initial value problems. We then 
show the Vaidya spacetime to be a limiting case of the 
general fluid equations, and the Schwarzschild spacetime 
to be a further limiting case of this. 

Due to the method of formulation, we can prescribe 
data on an initial spacelike hypersurface, apply physi- 
cally reasonable conditions, and the evolution of the sys- 
tem will take care of itself. In particular, we discuss at 
length a scenario whereby the initial hypersurface is sim- 
ply a ball of fluid extending out to some finite radius, 
and then vacuum extending beyond this to spacelike in- 
finity. As the system evolves, null radiation is allowed to 
emanate from the collapsing ball of fluid. Thus, at any 
spacelike slice beyond the initial hypersurface there will 
be three regions of spacetime; a collapsing fluid ball, a 
region of outgoing null radiation and a Schwarzschild ex- 
terior extending to spacelike infinity. Furthermore, as the 
Vaidya and Schwarzschild regions are shown to be limit- 
ing cases of the general fluid ball, the entire spacetime is 
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described in terms of a single line element, expressed in 
a single coordinate patch. 

As an expanding model is simply the time reverse of 
a collapsing scenario, the above model also satisfies ex- 
panding cases. The limiting case of such collapsing or 
expanding models however is a static model. As a final 
application we show the conditions for the set of equa- 
tions derived herein to be static. Using this we show a 
static spacetime necessarily has zero "heat flux" , and fur- 
thermore show the remaining equations to generalize the 
Tolman-Oppenheimer-Volkoff (TOV) equations of hydro- 
static equilibrium to include anisotropic stresses [H, [13] • 

The article is set out as follows; in section[H]we review 
the necessary components of the 3 + 1 formalism, and de- 
rive the equations governing the general fluid system. In 
section Uni we show how the equations are generalizations 
of both the Vaidya and Schwarzschild spacetimes, and 
pick out specific examples of each. Section IIVI discusses 
the overall picture of the gravitational collapse of a ball 
of fluid which emanates null radiation. Finally, in section 
W\ we show the reduction of the general equations to the 
static system. Throughout the article we use a +2 signa- 
ture and coordinates are denoted 2+ = (t, r, 6, </>) unless 
otherwise stated. Geometrized units are used, whereby 
c = G = 1. Greek indices range from ... 3, Latin from 
1 ... 3 and all other conventions follow @ ■ 

II. GENERAL FLUID DERIVATION 

We provide here a derivation of the general spheri- 
cally symmetric fluid equations, which closely follows the 
method used in [12j. To begin, we define a hypersurface 
forming, normalized, timelike vector field, n M , which has 
components 

n"=i(-l,/3,0,0). (1) 
a 

Here, a(t,r) > is the lapse function and (3(t,r) is 
the radial component of the shift vector. The energy- 
momentum tensor can then be irreducibly decomposed 
into components on and orthogonal to the spacelike hy- 
persurface formed by the normal vector. This decompo- 
sition is given by 

T„ v = {p + P) + Pg^v + 2 J(m n v)+ !W ( 2 ) 

Here, p(t, r) and P(t, r) are respectively the energy den- 
sity and isotropic pressure, 

3n = \Pj{t,r),j{t,r), 0, 0] , (3) 

is associated with the "heat flux" and Ii^ u is the trace- 
free anisotropic stress tensor. We note here that one can 
further discuss each of these terms in more detail with re- 
gards to properties of the fluids [see the brief discussion 
regarding equation (|28p below], however, the important 
aspect here is that this is a unique, irreducible decompo- 
sition of the energy momentum tensor. 



Following [12j , we define a mixed tensor according to 

Pi j :=diag[-2,l,l]. (4) 

Now, consider some (2) -tensor, say, which is sym- 
metric, spatial {W fjLa n a = 0) and trace-free (W a a = 0). 
It is trivial to show with spherical symmetry we can al- 
ways write 

Wy=tu(t,r)i^, (5) 

where w(t,r) is the distinct eigenvalue of Wij. With 
spherical symmetry, many of the 3 + 1 variables can be 
expressed in terms of their distinct eigenvalues, and we 
therefore make the following definitions; 

• The anisotropic stress tensor 

II,,: Hf'.ri/',,. (6) 

• The trace-free extrinsic curvature 

K ij --± ij K:=a(t,r)P ij . (7) 

• The trace-free three-Riemann tensor 

*R ij --L ij z R:=q{t,r)P ij . (8) 

• The trace-free Hessian of the lapse function 

~ \ DiDj - ~ +y D k D k )j a := e(t, r)P ir (9) 

• The electric curvature tensor 

Eij := C lajf3 n a n p := X(t, r)P l3 . (10) 

Here, and K are the extrinsic curvature and its trace 
of the spatial three-slices respectively, +y is the three- 
metric on the three slices, 3 Rij and 3 R are the three-Ricci 
tensor and scalar respectively, Di is the unique three- 
metric connection and C^^r is the Weyl conformal cur- 
vature tensor. 

As a final generalization to we must include one 
more non-trivial metric variable in the line element. We 
do this as the energy momentum tensor now has four 
degrees of freedom (i.e. p, P, j and n), implying to keep 
the fluid as general as possible we require a corresponding 
number of degrees of freedom in the metric. Therefore, 
we designate the spherically symmetric line element to 
have the general form 

dS * = .jo* + m + drf R2dQ2 

1 + E 

Here, E(t,r) > —1 is an otherwise arbitrary function 
which reduces to the energy function in the LTB dust 
model [ll[ and dfl 2 := d6 2 + sin 2 8d(j> 2 . Furthermore, 
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R = R(t, r) is an arbitrary function whereby, without 
loss of generality, we set R(t, 0) = 0. 

The EFEs can now be expressed in 3 + 1 form, and 
in spherical symmetry are given by two non-trivial con- 
straint and two non-trivial evolution equations. The 
Hamiltonian and momentum constraints are respectively 



? P + -K 2 -6a 2 =167rp, 
o 



Or 



3 



-K 



3a dR 
p" d7 



= - 4:1V j, 



(12) 
(13) 



and the two evolution equations are 



2C n K - - 3 R - K 2 - 9a 2 + -D k D k a =24ttP, (14) 
2 a 



C n a — aK + e — q = 87rII. 



(15) 



Here, C n is the Lie derivative operator with respect to 
the normal vector which, when acting on a scalar takes 
the form 



ld_ pd_ 

a dt a dr 



(16) 



By putting the line element through the momentum 
constraint (1131). one can show 



Regularity of p and j at r — 0, and given R(t,0) = 
implies M(t, 0) = 0. We refer to this function as a 
mass function as in the dust case it reduces to the mass 
function of the LTB metric [ll[ , and to the Schwarzschild 
mass in the vacuum case. Furthermore, we will show as 
one of our main results that it also reduces to the mass 
function of the Vaidya system (see section [ill B|) . 
Now, integrating both sides of equation (fT9|) gives 



M , , 1 + EdadR 

+ 4tt P - 211 )R = — — 

R z a or or 



(21) 



where the specification of R(t, 0) = and regularity at 
r = ensures any functions of integration vanish. 

Now, putting equation (|2"Tj) back through (fig)) , we find 
an equation for the square of C n R, which is 



(CnRY 



2M 
~R 



+ (l + E) 



dR 
dr 



(22) 



We now substitute this into equation (f2"Tj) . and also uti- 
lize equation (|17p to find the Lie derivative for the mass 
function 



C n M = -AirR 2 



(P-2U) C n R + j^(l+E) 



(23) 



OR 
dr 



2(1 + E) a dr 



= AivjR-^{L n R). (17) 



To deal with the evolution equations, we take equation 
(fT4|) and subtract six times equation (fl~5|) . That is, we 
obtain an equation giving the Lie derivative of the vari- 
able a + K/3. Expressing this equation in terms of the 
metric coefficients gives, after some algebra, 



1-(1 + E) 



dR 



dr 

2R(1+E) da dR 
a dr dr 



8tt (P - 211) R 2 



2R£ n 2 R-(£ n R)\ (18) 



where C n 2 R :— C n (C n R). It is both interesting and per- 
tinent to note that the two evolution equations, (|14[) and 
(|15[) . are dependant on one another due to the spherical 
symmetry. That is, these two equations can be reduced, 
without loss of generality, to the single evolution equation 
given by (fTS]) . Substituting this through the Hamiltonian 
constraint (|12j) . and utilizing the momentum constraint 
(fTB")) . gives the following equation, 



d_ 

dr 



R 2 (1 + E) da dR 



47T p 



a dr dr 
dR 
dr 



R 2 C n 2 R - 4tt (P - 2n) P 3 
]C n R ) R 2 . (19) 



We now define a function, M(t,r), according to 

dM f dR „ \ 2 

^r:=4J P - +J C n R)R 2 . 



(20) 



The conservation of energy- momentum equations, 
VaT"^ = 0, provide the remaining pieces of information 
regarding the spacctimc. When decomposed onto and 
orthogonal to the spacelike hypersurfaces these become 

= (Cn -K)p + D k j k + 2h k Jk -PK- ma, (24) 
= (C n - K)ji + (D k + n k ) (U lk + L lk P) + phi, (25) 

where n, is the four-acceleration. By putting the line ele- 
ment through these equations, they respectively become 



C nP - (p + P - 2U) 
'Idj , 



C n E 



J (1 + E) 
—C n R, 



2(1 + E) 
1 dE 



j dr 2(1 + E) dr 



1 d(3 2 

a dr R 

2 dR 2 da 

R dr a dr 



C n R 



(26) 



C n (jy/T+E) _ 2 , 

VTTe 3 

(p + P - 211) da 

a dr 



1 



2(1 + E) " + adr R n 

o .„ 6nap 



The two conservation equations expressed above provide 
evolution equations for two of the four matter variables, 
namely p and j. Therefore, the only way of evolving 
the other two matter variables, P and n, is to employ 
two equations of state to supplement the system. Fur- 
thermore, we note the system of equations is currently 
underdetermined. One simple way of expressing this is 
to note the differences between the general fluid system 
and the more restrictive perfect fluid system presented 
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in [l2|, which was fully determined. The key differences 
between the two systems are that we now have two more 
matter variables, j and II, and one more geometric vari- 
ables, R, however, only two extra equations have been 
introduced, which are (|26[> 1 and one more equation of 
state. 

The system being underdetermined implies that to find 
solutions for various physical scenarios requires the input 
of extra physical or geometric effects. One possible way 
is to introduce an equation governing the conservation 
of Baryon number. In this way, Baryon number will be 
conserved for the interior fluid region, implying only non- 
baryonic energy may escape the collapse region. This 
conservation equation is usually expressed as (see for ex- 
ample [IE 01) 



V Q { Po u a ) = 0, 



(28) 



where po is the rest-mass energy density, and u M is the 
comoving four-vector. The introduction of this equation 
will close the system, however its introduction requires 
the additional consideration and discussion of thermody- 
namic quantities, a task reserved for future work. 

Summarily, the system of equations describing a gen- 
eral spherically symmetric fluid are the line clement pip 
and the equations governing the metric coefficients (|17[) . 
(f22|) and (|23| , along with two equations of state and the 
two conservation equations, (f2l>)) and (f2"T|) . 



A. Tidal Forces 

It is an interesting exercise to calculate the tidal forces 
for the general fluid spacetime, and thereafter for the 
resulting subclasses of Vaidya and Schwarzschild space- 
times. The tidal forces are given by the electric curvature 
tensor [13] , E^j , which is given by the contraction of two 
indices of the Weyl tensor with respect to the normal 
vector [see equation (1101) ]. A constraint equation from 
the gravito-electromagnetic system of equations, when 
expressed in terms of metric coefficients, is 2 



9 A tt 4 ^ 



3 dR 
R~dr~ 



(A + 4nU) 



4ttj 
R 



C n R. 

(29) 



Substituting equation (|20p into the above, we find an 
equation which can be integrated to give 



A 



M 
R? 



4tt 



(p + 3n). 



(30) 



Here, we have used R(t, 0) = and regularity of A, p 
and II at r = 0. We will see in the following sections 



1 In the perfect fluid case, this equation is an identity which is the 
derivative of equation (123 D - 

2 The magnetic component of the Weyl tensor vanishes identically 
in spherical symmetry. 



how this becomes even simpler for both the Vaidya and 
Schwarzschild spacetimes. 



III. EXTERIOR REGIONS 

Having derived the fluid equations we can show how, 
by specifying particular equations of state, these reduce 
to both the Schwarzschild and Vaidya spacetimes. This 
implies we can establish collapse scenarios by having 
an interior fluid, and an exterior region being cither 
Schwarzschild or Vaidya. Moreover, the entire spacetime 
will be described in a single coordinate system, by simply 
taking the interior metric and letting the matter variables 
tend to the appropriate equation of state at some finite 
radius. 



A. Schwarzschild 

The reduction to the Schwarzschild spacetime is simply 
given by having all the energy-momentum variables van- 
ish at some radius (i.e. p = P = j = H = 0). Equation 
([2H)l and (|2"3")l then imply that the mass function is a con- 
stant, i.e. M — M s . The tidal forces for the spacetime 
are now given by 



A = 



i? 3 



(31) 



This is a slightly more general result than the standard 
form (see for example [llj) due to the function R deter- 
mining the sizes of the two-spheres. 

There now remains four metric coefficients, a, /?, E 
and R, and just two non-trivial equations, which are 



dR 
dr 



C n E 



2(1 



(C n RY 



E) 
2M B 
R 



180 
a dr 



= --(C n R), (32) 



Therefore, the line element (fTT| where a, 0, E and R 
satisfy equations (|3"2"]) and (f3"3"|) describe an infinite fam- 
ily of coordinates for the Schwarzschild spacetime. This 
system is underdetermined as there are two differential 
equations for four variables, and as such, one may make 
any suitable choice for any of the four metric coefficients, 
provided equations (|32|) and ((33)) can be satisfied. 

Whilst this is a perversely general system describing 
the relatively simple Schwazrschild spacetime, we should 
note that the amount of generality is a necessary evil. 
This is because, in order to describe the entire system in 
a single coordinate patch as per our initial aim, one must 
first establish the interior, matter filled region, and then 
use the coordinate freedom in the exterior to express the 
Schwarzschild region in appropriate coordinates. Thus, 
the amount of freedom in the exterior corresponds to the 
vast number of solutions for the general fluid equations 
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for the interior. As an example, if we consider matching 
a general fluid to a vacuum region (with no intermediate 
Vaidya region) then, as in the case of [l2j, we demand 
there exist a small region of dust between the fluid and 
the vacuum. This dust region has all matter terms van- 
ish with the exception of p. Equation (|27|> then implies 
the lapse is only a function of the temporal coordinate, 
and hence coordinate freedom can be utilized to set it to 
unity without loss of generality. Therefore, to match the 
vacuum region to the dust region, we also set a = 1 in 
the Schwarzschild region. 

With a = 1 for the Schwarzschild region, many familiar 
solutions arise. This coordinate choice corresponds to a 
vanishing four-acceleration, n l = 0, and thus the normal 
vector is tangent to a congruence of timelike geodesies. 
If we make the further choice R = r, then the system 
reduces to the familiar generalized Painleve-Gullstrand 
(GPG) coordinates [THEl]. Another possibility is to set 
(3 = 0, which implies equation (f3"2)) can be integrated and 
we find the Novikov coordinates 



dS 2 



-dt 2 



dR 
dr 



dr 2 + R 2 dfl 2 



(34) 



We have denoted R to be the radial coordinate of the 
Vaidya metric such that we can perform a coordinate 
transformation given by R = R(t,r) and u = u(t, r), 
where t and r are the temporal and radial coordinates of 
the general fluid line element respectively. We define the 
simultaneous coordinate transformation to be 



dR _ am (3(l-m/R) dR_-{l-m/R) 
~dt ~ R VTTe ' ~dT ~ y/TTE 
du (3 du 1 

~dt = ~ a + TJTTE z, = VTTe' 



(38) 



which renders the line element in the exact form of equa- 
tion (fTT|) . The above equations imply the important re- 
lation 



(39) 



Furthermore, the integrability conditions must be sat- 
isfied for this to be a valid coordinate transformation. 
These integrability conditions provide restrictions for the 
metric coefficients which, after some algebra, are found 
to be 



where / = f(r) and R is a solution of 

2 



OR 
dt 



2M S 
R 



(35) 



We have shown how a fluid matches directly to a 
Schwarzschild region. However, our aim was to include 
radiation into the system, which implies we will actually 
be matching the fluid to a Vaidya region. Furthermore, 
the Vaidya region will match onto a Schwarzschild region 
along a null line. To establish these conditions, we must 
first establish the reduction of the general fluid equations 
to the Vaidya system. 



B. Vaidya 

It is not intuitively obvious how one retrieves the 
Vaidya spacetime from the system of equations derived 
herein. This is because we do not know a priori appro- 
priate equations of state for the Vaidya spacetime when 
written using non-null coordinates. For this reason we 
derive the results of this section in reverse. That is, we 
start with the Vaidya line element in standard coordi- 
nates, (u, R, 9, 4>), which is given by 



dS 2 



2m 
~~R 



du 2 



2dudR + R 2 dfl 2 , (36) 



where u is an outgoing null coordinate and m{u) is the 
Vaidya mass function satisfying 



dm o 
— = AirpR 2 
du 



(37) 



VTTE da C n E 



a 



dr 



2(1 

-1 



-E) 
C n m 



ld£ 
a dr 

-VTTe 



dm 
dr 



m 
R 2 ' 



(40) 



Therefore, the Vaidya spacetime is given by the line 
element ([TTT) . whereby the metric coefficients satisfy 
It is straightforward to see this 
as with the 



equations (|38M0p 

system is underdetermined and therefore, 
Schwarzschild case, there are infinitely many ways of ex- 
pressing the Vaidya spacetime in these general spheri- 
cally symmetric coordinates. That is, any solution of 
the above equations represents a Vaidya spacetime, the 
choice of which depends on the physical problem being 
studied. 

We can now put these conditions through the reduced 
field equations we have derived to obtain the equation of 
state for the Vaidya spacetime when expressed in non- 
null coordinates. Beginning with equation (|22p we find 
the remarkable result that 



m = M. 



(41) 



That is, the mass we derived with the definition given by 
equation (|20|) reduces to the Vaidya mass in the appropri- 
ate limit. This function was also shown in [121 ] to reduce 
to the Schwarzschild and LTB masses in the appropriate 
limits. 

Now, if we put the integrability conditions through 
equations (jTTJ) and PTj) . after some algebra we can show 



C n M + VTT~E^- = AttR 2 (p - 211 - jVTT~Ej 



(42) 
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Furthermore, putting the integrability conditions 
through equations ([20")) and ([25]) , we find respectively 



dM _ -AttR 2 
~dr ~ VTTE 



M / , 

-^{p-jVTTe 



(43) 



C n M =AttR 2 



— (p-2U- jVT+E) + j \ I E 



(44) 



Finally, combining equations (|42j|44p , we derive an alge- 
braic relation between all energy momentum variables 



P- 



P - 2H - 2jVl + E = 0. 



(45) 



Now, putting the above information into the conservation 
equations (|26|27|) . and combining these we find 



c n + VT+e^-) (p - jVT+E 



R 



iVT+E + 3n 



(46) 



The remaining piece of information to be utilized is the 
form of the mass in the Vaidya solution, i.e. equation 
(|37|) . By inverting the coordinate transformation and 
expanding the derivative with respect to the null coordi- 
nate in terms of radial and temporal derivatives, one can 
show 



dm 
du 



4nR 2 jVT+E, 



(47) 



where we have also made use of equations (|42M5|) . 
Hence, direct comparison with equation (|37p implies 
p = j\fl + E. Finally, combining all information, we 
find the equations of state governing the Vaidya space- 
time are 



Also putting the equations of state back through equation 
(HOI , we find the lapse function is highly constrained 



1 da 
a dr 



-M 



R 2 VT+E 



(51) 



This equation has some interesting physics in that the 
right hand side only vanishes if the mass function van- 
ishes. Therefore, providing M ^ 0, the lapse function 
must necessarily be a function of the radial coordinate, 
and hence we are never free to set the lapse function to 
unity for the Vaidya spacetime. Another way to state this 
result is that in the Vaidya spacetime, the unit timelike 
vector normal to spacelike slices, n M , can not be tangent 
to timelike geodesies. 

Furthermore, this equation gives a matching condi- 
tion for the Vaidya and Schwarzschild regions. As men- 
tioned, the Schwarzschild system is vastly underdeter- 
mined, however, when matching Schwarzschild to Vaidya 
we require the lapse be a solution of (|5ip in the vacuum 
region (with M = M s ) such that it matches continuously 
across the interface. This point is discussed in more de- 
tail in section HV CI 

As mentioned, the Vaidya system of equations is un- 
derdetermined. It is interesting to note that the choice 
of functions R ~ r is not a valid one for the Vaidya 
spacetime. This is because choosing R — r reduces the 
system to two differential equations, the general solution 
of which has the mass function being a constant. Hence, 
choosing R — r forces the Vaidya system to give the 
Schwarzschild solution. Therefore, establishing a radi- 
ating problem necessitates the extra degree of freedom 
we put into the line element in the form of the general 
function R. 



3P = jVT+E 



-3n. 



(48) 



It is interesting that p = 3P is a radiative equation of 
state, and hence this aspect of the result is one that is 
expected. 

We can put this equation of state through equation 
([30| to show the tidal forces for the Vaidya region of the 
spacetime are given by the simple result 



(49) 



This is the same form as for the Schwarzschild spacetime, 
however we note the mass in the Schwarzschild is a con- 
stant, whereas it is still a function of both the temporal 
and radial coordinates in the Vaidya case. 

There are two further results from this section that 
will prove fruitful. Firstly, putting the equations of state 
through equation (|42|) we find 



-dM 



C n M + VT+E— = 0. 

or 



(50) 



IV. COMPLETE PICTURE 

A. Summary of Equations 

We provide here a summary of the important equations 
for the general fluid system, as well as for the reduction 
to both the Vaidya and Schwarzschild spacetimes. This 
enables us to see the appropriate limits one needs to take 
at the boundaries between any two regions. 



General Fluid 



There are three main evolution equations resulting 
from EFEs. These, and the mass equation, are given 
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by 



R'C n E R'd/3 d 

WTm = AnjR ~ flF " 5F ' (52a) 

(£„i?) 2 =^ + (l + £)i?' 2 -l, (52b) 



£„M = - AttR 2 [(P - 2n) C n R + jR' (1 + E)] , 
<9M 



c9r 



=4tt (pR' + jC n R) R 2 , 



(52c) 
(52d) 



where R' := dR/dr. The two conservation equations, 
(J2BJ) and l[2Tp. are also required for the system. 



2. Vaidya 

The Vaidya system is given by the equation of state 
(|48[). and the same three evolution equations and mass 
definition become 



R'C n E _ AvpR R'dp d 
2{l + E) ^TTE a 8r Or ( n ' 

(C n R) 2 ~- + (1 + E) R' 2 -l, 
K 

C n M =4npR 2 , 
dM —AnpR 2 
~~dr ~ y/T+E" 



(53a) 

(53b) 
(53c) 

(53d) 



Now, the two conservation equations have reduced to a 
single equation, 



C n p + VI + E— = — 
or R 



particular aspects of the entire system [including equa- 
tions (f5"3"|) and ([54")) ]. In particular, we consider a sce- 
nario whereby we define the data on an initial spacelike 
hypersurface (denoted So in figure [T]) such that a ball of 
spherically symmetric fluid reaches out to a finite radius 
of r = r*(i = 0), such that vacuum exists for r > r*(0). 

We then allow the system to evolve, with the condition 
that the interior fluid ball emanates only null radiation. 
At present this condition is put in manually by asserting 
the boundary of the fluid region to be where the equations 
of state for the interior reduce to the equation of state 
for the Vaidya system, i.e. equation (|4"5|) . We therefore 
have the boundary, r*(t), separating the region of general 
fluid, region I of figure[U and the Vaidya region, denoted 
region III. This Vaidya region will have begun on the 
initial hypersurface with zero volume, and be expanding 
along the null cone in the direction of increasing radius, 
extending out to future null infinity, denoted J zr+ . Ex- 
terior to this Vaidya region is a region of Schwarzschild 
spacetime, denoted as region IV and governed by equa- 
tions (|54p . extending out to spacelike infinity, i. 



(53e) 




3. Schwarzschild 

For the Schwarzschild region, all the matter terms van- 
ish, and hence the evolution equations are 



R'C n E _ R_dl3_ _ d_ 
2{l + E) ~^~dr~ ~ ~b~r~ 



(C n R) , (54a) 



(C n Rf =^f + (l + E)R' 2 -1, (54b) 



C n M =0, 
dM 



dr 



0, 



(54c) 
(54d) 



where the conservation equations are now trivially satis- 
fied. 



B. Generic Collapse Model 

While analytic solutions of the general fluid equations 
are obviously difficult to find, we can still analyse 



FIG. 1: Compactification diagram for the collapse of a ball of 
fluid: i represents spacelike infinity, is future null infinity, 
I + is future timelike infinity and Eo is the initial spacelike 
hypersurface. 



We note the other condition we have imposed for this 
scenario is that we used an equation of state and initial 
conditions such that the radius of the fluid region reduces 
to zero volume in a finite time. Using equation (|52H54[) . 



we can look at what happens when this collapse point is 
reached. In particular, when r* — > 0, we note the energy 
density will diverge. The mass function will then become 
a constant, on the proviso that the j does not diverge. 
The mass function becoming a constant is enough to im- 
ply the system of general fluid equations (|52|) reduce to 
the system of Schwarzschild equations (|54[) . Therefore, 
as soon as r* = a Schwarzschild black hole remains, as 
indicated by figure Q] 

We further look at the apparent horizon, which is indi- 
cated schematically by the dashed line in figure [T] This 
will evolve from r = at a finite time, and eventually 
intersect with the surface of the fluid ball, r = r*. The 



8 



interior fluid region will continue from this point unaf- 
fected, however, the exterior will no longer be the Vaidya 
system. This is due to the surface, r*, being inside the 
apparent horizon, implying the null cones have tilted over 
to such an extent that null radiation can no longer escape. 
Therefore, the region interior of the event horizon, but 
still exterior to the surface r = r* is Schwarzschild (re- 
gion II) . This event horizon extends out to future timelike 
infinity, X + . 

As a final aid, figure [2] shows a schematic of spacelike 
slices through the compactification diagram. We note 
that each slice, schematically represented by the dashed 
lines, begins in the general fluid region, extends through 
Vaidya, into the Schwarzschild region and out to space- 
like infinity, i. 




FIG. 2: Compactification diagram showing the evolution of 
spacelike hypersurfaces 



C. Matching example 

Specific solutions of the model we have established will 
usually be found by first constructing the fluid region, 
and using this to determine the form of the regions ex- 
terior to it, i.e. the Vaidya and Schwarzschild regions. 
However, analytic solutions for the interior are difficult 
to come by and will probably require numerical analy- 
sis. Therefore, as an illustrative example, we find here a 
particular example for the Vaidya region, and show how 
this reduces to the Schwarzschild region in the same co- 
ordinates. In particular, the example we choose is the 
marginally bound case, such that E — 0. This is the re- 
verse scenario to how one would tackle a physical prob- 
lem, and we shall see that this implies the Vaidya region 
is underdetermined, as it requires the input of boundary 
conditions coming from the fluid region. 

The amount of coordinate freedom in the vacuum re- 
gion implies we can stipulate conditions on the metric 
functions according to the form of these functions in 
the Vaidya region. That is, in the Vaidya region, equa- 
tions (|3"8")1 and (|5ip (now with E = 0) hold, so without 



loss of generality we can stipulate these to hold for the 
Schwarzschild region also. This will imply both regions 
will be expressed in a single coordinate patch, as per our 
initial aim. In the Schwarzschild region the mass becomes 
a constant, and equations (f38|) and (|5Tj) become 



dR 
dr 
dR 
~dt 
1 da 
a dr 



1 



= -/m- 

M. 
R 2 ' 



M, 

M, 
~R 



aM s 
R 



(55) 
(56) 
(57) 



Furthermore, simplifying equation (|54a| with E = and 
the above equations gives 



180 
a dr 



R 2 



(58) 



These four equations can be solved to give an exact solu- 
tion for the exterior Schwarzschild region when matched 
on to a marginally bound Vaidya region. Firstly, to inte- 
grate equation (|55[) we require the Lambert W function, 
W(£), which is defined such that it is a solution of the 
algebraic equation (see for example [l8l |) 



Then, by defining 



■/(*)] 



(59) 



(60) 



where f(t) is an arbitrary function of integration, we find 
the solution of equation (|55p is 



R = M s [W(0 + 1] 



(61) 



The function f(t) is given by utilizing a more spe- 
cific boundary condition between the Vaidya and 
Schwarzschild regions. It is interesting to note at the 
end point of collapse (when the entire solution is sim- 
ply given by the Schwarzschild spacetime), we can use 
the boundary condition that R(t, 0) = 0, which implies 
f(t) = — M s lnM s , which further simplifies the form of 

Now, putting (|6ip through the remaining three equa- 
tions, we find 



= 



-i df 

W{0 + 1 + di' 



(62) 
(63) 



where we have used coordinate freedom to re-scale an 
arbitrary function of t to unity in the solution of the 
lapse function. Thus, the line element (fTTj) , where the 
metric coefficients are given by equations (I6H63P with 
E = 0, describes the Schwarzschild spacetime in a coordi- 
nate system which matches continuously to a marginally 
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bound, non-null coordinate system for the Vaidya space- 
time. 

While we have solved for the exterior Schwarzschild 
region, this does not provide enough information to solve 
for the Vaidya region. We can show in the Vaidya region 
that equations (|5Tj) and (|53bj) imply 



a — [3 — 



df 



(64) 



which further implies equation (|50|) can be integrated to 
give 



M = M [r + f(t)} . 



(65) 



If we integrate equations ([55]) for the null coordinate u, 
we find u = r + /(t), which implies the expected re- 
sult that the mass is a function of the null coordinate, 
M = M(u). Therefore, characteristic curves of the mass 
function are the null lines in the spacetime. Moreover, 
this implies that the boundary between the Vaidya and 
Schwarzschild region, ro(£), is given by a specific charac- 
teristic, implying this boundary is null. 

In order to integrate the remainder of the system en- 
tirely, we need to solve 



dR_ ( M 
7fr ~ ~ { ~~R 



(66) 



which is an Abel equation of the second kind. However, 
to integrate it we require more knowledge of the spe- 
cific form of the function M, which is to be gained from 
boundary conditions. Unfortunately, the boundary at 
ro(t) offers no new information on the form of the mass 
function. Intuitively, this is due to the causality of the 
system, which essentially works from the interior to the 
exterior regions (i.e. the interior regions influence the 
exterior regions). Therefore, the boundary conditions re- 
quired to solve the remainder of the system are given by 
the boundary at r*. Hence, once a specific fluid region 
is found, it can be used to stipulate boundary conditions 
at rv, which can then be used to fully solve equation 
([6"6"| . Once this equation is solved, equation (|5"Tj) can be 
integrated to find the lapse function, equation (|64[) im- 
plies we know /3, and we can use the remaining equation 
on dR/dt to relate the arbitrary functions. Therefore, 
the knowledge of the functional form of M, given by the 
boundary conditions at r*, will imply the Vaidya region 
is completely soluble. 



V. STATIC SOLUTION 

There is the situation whereby a model is not collaps- 
ing or expanding, but is static. In reality, while a star 
may appear static, it will still radiate, and hence is not 
static in the relativistic sense. A purely static system 
however will have zero fluxes, implying it will not radiate 
and hence the exterior spacetime is simply Schwarzschild. 
The compactification diagram for this scenario is given 



in figure [3l and this section is therefore devoted to using 
our derivation of the general fluid equations to find con- 
ditions for the spacetime to be static. In the process we 
generalize the TOV equations of hydrostatic equilibrium. 
We note however that there should further exist solutions 
to our equations that are not static, but the boundary of 
the fluid region remains constant. These solutions will 
be sought in a future article. 




FIG. 3: Compactification diagram for static scenarios 

A static spacetime is defined such that there exists a 
timelike, Killing vector field which is everywhere hyper- 
surface orthogonal . The normal vector we have used 
is both timelike and hypersurface orthogonal, however 
it being normalized places too tight a restriction on the 
function. We therefore generalize the normal vector by 
multiplying by an arbitrary function, such that 



(67) 



where T = J-(t,r). It is trivial to show is also 
timelike and hypersurface forming. Therefore, solving 
Killings equations with respect to this vector field will 
provide conditions for staticity of the spacetime. 
Killings equations with respect to iV M are 



dx a 



8N a 
dx u 



9au 



dN a 
dxf* 



0. 



(68) 



which provides a system of ten coupled differential equa- 
tions. While a majority are trivially satisfied, we find 
evaluating the (6, 9) or (0, cf>) components gives 



C n R = 0. 

Likewise, the (r, r) component implies 

CnE , 18(3 



2(1 + B) 



dr 



= 0. 



(69) 



(70) 



Evaluating the (t, r) component, and utilizing equations 
(|69|) and ([70]) gives a condition on T with respect to the 
lapse function 



d_ 

dr 



= 0. 



(71) 



10 



Integrating this gives a = g(t)J- ', where g(t) is an arbi- 
trary function of integration. One can always utilize co- 
ordinate freedom to rescale the temporal coordinate such 
that g(t) scales to unity. Therefore, we find the relation 
T = a, implying the Killing vector is given by 



d_ 

dt 



d_ 
dr' 



(72) 



Finally utilizing equation (|69[) along with the (t, r) com- 
ponent of Killings equations gives 



C n a = 0. 



(73) 



Therefore, if the metric coefficients satisfy equations 
(|69|) . ([70]) and ([73]) . then the spacetime admits a hyper- 
surface orthogonal, timelike Killing vector field, imply- 
ing the spacetime is static. We note that the standard 
method for searching for static solutions in a spacetime 
is to express the line element in a diagonal form, however 
we see the conditions given by the equations ([S"5]) . ([70]) 
and (|75|) generally give (3 ^ 0, implying the static space- 
time is written in non-diagonal form. One can further see 
this implies d/dt is not the Killing vector, and in general 
there will be dependance on the temporal coordinate in 
the line element. 

Putting ([nU) and ([73"]) into equation (fT7|) instantly gives 
the condition 



i = o. 



(74) 



That is, a static spherically symmetric fluid necessarily 
has vanishing "heat flux" . This makes physical sense as 
a flux intuitively implies a time rate of change, which is 
ruled out by the staticity condition. 
Now, equation ([2"2"]1 reduces to 



and furthermore from (123|) we find 

C n M = 0, 

where the mass function is now 
dM . dR 



2M 



4tt P ^R 2 . 
or 



(75) 



(76) 



(77) 



Putting these results through equation ([2"Tj) and l[2"7|) 
gives the generalization of the TOV equation for hydro- 
static equilibrium 



1 da 

a dr 



M + An (P - 211) R 3 dR 



R 2 (1 - 2M/R) 
-1 



dr 



P-2U 



d , . 6n dR 



(78) 



These equations are a generalization of those found in fl3l . 
[l4j], whereby we have extra freedom in the coordinates 



in form of the functions (3 and R. Finally, putting the 
staticity conditions through equation (|2"6"|) gives 



C n p = 0. 



(79) 



Summarily, the static system is given by equations 
HS1]), (HQ]), C3l and ^MB- 
it is trivial to show if we let R = r, then equation 
([55]) implies (3 = and the line element reduces to the 
standard diagonal form. One can then see that a, E, p 
and M all become functions of only the radial coordinate, 
which is consistent with d/dt being the Killing vector. 
Furthermore, the equations directly reduce to those given 
in [3, 03] , who also provide some exact solutions of the 
equations by specifying various equations of state (also 
see [13] and references therein). 



A. Exterior static region 

As with the non-static case, we can utilize the form of 
the interior to find the form of the exterior Schwarzschild 
region which matches continuously to the interior. We 
begin by letting all the matter variables vanish at some 
finite radius, which implies M = M s . The left hand equa- 
tion in ([78]) can then be integrated for the lapse function 
to give 



a = \ 1 - 



2M S 
R 



(80) 



where coordinate freedom has again been used to set a 
function of integration to unity. It is interesting that 
there still exists coordinate freedom for the exterior, 
which would be determined when matching to a more 
specific form of the interior. We choose here, as an illus- 
trative example, to let R = r. As mentioned, this implies 
[3 = 0, and equations ([70]) and ([75]) can then be used to 
show 



E= — 



2M S 



(81) 



This implies the line element becomes the familiar 
Schwarzschild metric 



dS 2 



2M, 



dt 2 



dr 2 



1 - 2M s /r 



r 2 dn 2 . (82) 



Therefore, a static interior general fluid, with R = r, 
matches continuously onto a Schwarzschild region in 
Schwarzschild coordinates. 



VI. CONCLUSION 

We have derived general equations governing a spher- 
ically symmetric system with an arbitrary matter dis- 
tribution. We showed that the system is underdeter- 
mined, and requires the input of more physics, most 
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likely in the form of an extra conservation equation and 
thermodynamic considerations. Furthermore, we have 
shown how the system reduces both to the Vaidya and 
Schwarzschild spacetimes by imposing only particular 
equations of state. This implies entire collapse problems 
can be established utilizing the same coordinate patch 
throughout the spacetime. This provides an alternative 
approach to the study of gravitational collapse of gen- 
eral fluids, which is usually achieved by establishing an 
interior metric that joins to the Schwarzschild or Vaidya 
metrics in the exterior usingiunction conditions at the in- 
terface (see for example [1, ll, @] ) ■ However, this assumes 
one knows the form of the exterior metric a priori. With 
the inclusion of more physics into the method derived 
herein, one will be able to determine the exterior region 
of the spacetime without prior knowledge. 

As an example of this, we describe the equations for the 
set-up of an initial value problem whereby on an initial 
spacelike slice a ball of fluid extends out to some finite 
radius, beyond which is vacuum extending to spacelike in- 
finity. As the system evolves, the interior is permitted to 
emit null radiation, and hence a region of Vaidya space- 
time extends between the fluid and the vacuum. While 
this scenario was described qualitatively, specific models 
of these scenarios should be made possible with numerical 
solutions of the equations derived herein. This should al- 
low for the modelling of more physical problems, whereby 
one does not require prior knowledge of the final solution 
in order to solve the problem. 

This work opens avenues for further research into many 
currently unanswered questions. In particular, the intro- 
duction of thermodynamics into the system will allow for 



the inclusion of conservation laws, such that knowledge of 
the exterior region is not required a priori. Furthermore, 
being able to discuss the evolution of the temperature in 
all regions of the spacetime, including at r = 0, at the 
apparent horizon, and also at the boundary of the various 
regions of the spacetime is an exciting prospect. 

Another aim of the current program of research is to 
extend our method beyond that of spherical symmetry. 
In particular, we believe this method will be adaptable 
to the case of quasi-spherical symmetry, such that we 
can discuss interior regions for the Robinson- Trautman 
spacetime which represents outgoing gravitational radi- 
ation. We expect the compactification diagram for the 
quasi-spherical case will be of the same form as figure [TJ 
however region III will represent the Robinson- Trautman 
solution and region I will represent some form of interior 
fluid with deviations from spherical symmetry. 

Another future application of the current method is to 
introduce axisymmetry. It is still an open question as 
to whether the Kerr solution represents the exterior to a 
rotating system. Under the current paradigm, one would 
first derive general interior fluid equations, and simply 
let the matter terms tend to zero at some finite radius, 
thus determining the form of the exterior region of the 
spacetime. 
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